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ABSTRACT. We address the well-posedness theory for the magento-geostrophic equation, namely an active 
scalar equation in which the divergence-free drift velocity is one derivative more singular than the active scalar. 
In the presence of supercritical fractional diffusion given by (—A) 7 with < 7 < 1, we discover that for 
7 > 1/2 the equations are locally well-posed, while for 7 < 1/2 they are ill-posed, in the sense that there is no 
f~>) ■ Lipschitz solution map. The main reason for the striking loss of regularity when 7 goes below 1/2 is that the 

' constitutive law used to obtain the velocity from the active scalar is given by an unbounded Fourier multiplier 

which is both even and anisotropic. Lastly, we note that the anisotropy of the constitutive law for the velocity 
■ may be explored in order to obtain an improvement in the regularity of the solutions when the initial data and 

I the force have thin Fourier support, i.e. they are supported on a plane in frequency space. In particular, for such 

well-prepared data one may prove the local existence and uniqueness of solutions for all values of 7 £ (0, 1). 
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1. Introduction 

The geodynamo is the process by which the Earth's magnetic field is created and sustained by the mo- 
tion of the fluid core which is composed of a rapidly rotating, density stratified, electrically conducting 
fluid. Recently Moffat |[T9l proposed a model for the geodynamo which is a reduction of the full magne- 
tohydrodynamic system of PDE. This model can be viewed as a nonlinear active scalar equation in three 
dimensions. An explicit operator M [0] encodes the physics of the geodynamo and produces the divergence 
free drift velocity U(t, x) from the scalar "buoyancy" field <d(t, x). We call this active scalar equation the 
magnetogeostrophic equation (MG). It has some features in common with the much studied two dimen- 
sional surface quasigeostrophic equation (SQG) (see (H [3j 01 HI [141 [161 |2TJ [26| and references therein). 
However the (MG) equation has a number of novel and distinctive features due to the strongly singular and 
anisotropic nature of the operator M. 

The critically diffusive (MG) equation is defined by the following system 

d t e + u- ve = s + kAg (i.i) 

V • U = (1.2) 
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where k > is a physical parameter, and S(t, x) is a C°°-smooth, bounded in time source term. The 
velocity U is divergence-free, and is obtained from via 

Uj = Afj-0, (1.3) 

for all j £ {1, 2, 3}, where the Mj are Fourier multiplier operators with symbols given explicitly by 

O = 2^ 2 fe 3 |fc| 2 -(/3 2 /r ? )fc 1 fc 2 2 fc3 

~ -20fc 1 fe 3 |fc| 2 -(^/ ?7 )fclfc 3 

where the Fourier variable k G 7^ d is such that fc 3 / 0. Note that Mj(k) is not denned on k 3 = 0, since 
for the self-consistency of the model it is assumed that and U have zero vertical mean. It can be directly 
checked that kj ■ Mj(k) = for all k G Z d \ {/c 3 = 0}, and hence the velocity field U given by (11.3b is 
indeed divergence-free. The physical parameters of the geodynamo are the following: Q is the rotation rate 
of the Earth, r\ is the magnetic diffusivity of the fluid core, k is the thermal diffusivity, and (3 is the strength 
of the steady, uniform mean part of the magnetic field in the fluid core. The perturbation magnetic field 
vector b(t, x) is computed from 0(i, x) via the operator 

bj = (p/r))(-A)- 1 d 2 M j e, for all j e {1, 2, 3}. (1.7) 

We observe that the presence of the underlying magnetic field, reflected in the parameter /3 2 /r/, produces a 
non-isotropic structure in the symbols Mj. 

It is important to note that although the symbols M, are 0-order homogenous under the isotropic scaling 
k — > Xk, due to their anisotropy the symbols Mj are not bounded functions of k. To see this, note that 
whereas in the region of Fourier space where k\ < max{A:2, k%} the Mj are bounded by a constant, uni- 
formly in \k\, this is not the case on the "curved" frequency regions where k% = 0(1) and k 2 = 0(\ki\ r ), 
with < r < 1/2. In such regions the symbols are unbounded, since as | k\ \ — > oo we have 

|Mi(fci, M r , 1)| w |fci| r , |M 2 (fci, \h\ r , 1)| w \h\, \M 3 (h, \h\ r , 1)| w \h\ 2r . (1.8) 

In fact, it may be shown that |M"(fc)| < C\k\, where C(/3, rj, Q) > is a fixed constant, and this bound 
is sharp. The fact that symbols are at most linearly unbounded in the whole of Fourier space implies that 
M"[0] is the derivative of a singular integral operator acting on (see [9] for details), as opposed to the 
case of the (SQG) equation where U is the Riesz transform of 0. 

In ||9l we studied a class of nonlinear active scalar equations of which the (MG) system (ll.lb - dl.6b is a 
member, namely (11.1b with 

Uj = djTijG (1.9) 

for all j G {1, 2, 3} where Tjj is a 3 x 3 matrix of Calderon-Zygmund operators of convolution type such 
that didjTij(j) = for any smooth function <fi. Inspired by the proof of Caffarelli and Vasseur O for the 
global well-posedness of the critically diffusive (SQG) equations, we used DeGiorgi techniques to obtain 
the global well-posedness of the critically diffusive (MG) equations, namely (11-1 b — (1 1.6b - We remark that 
for both the diffusive (MG) and (SQG) equations criticality is defined with respect to the L°° maximum 
principle. In iflOl we considered the non-diffusive version of the (MG) equations, namely dl.lb - dl.6b with 
k = 0. In contrast with the critically diffusive problem where the (MG) equation is globally well-posed 
and the solutions are C°° smooth for positive time, we proved that when k = the equation is Hadamard 
ill-posed in any Sobolev spacd3 which embeds in W 1 ' 4 . 



^The non-diffusive (MG) equations are however locally well-posed in spaces of real-analytic functions, cf. 1101 . 
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In this present paper we study the fractionally diffusive version of the (MG) equation, namely 

d t e + U • VG = S - k(-A) 7 Q (1.10) 
V-U = 0, U = M[Q] (1.11) 

with Uj = Mj[@] where the symbols Mj are given by jl.4b - dl.6l > and the parameter 7 £ (0, 1). The main 
observation is that the structure of Mj (anisotropy, unboundedness, and evenness) and the weak smoothing 
effects of the nonlocal operator (—A) 7 combine to produce a sharp dichotomy across the value 7 = 1/2. 
More precisely, if 7 > 1/2 the equations are locally well-posed, while if 7 < 1/2 they are ill -posed in 
Sobolev spaces, in the sense of Hadamard. 

In section [3] we use energy estimates to prove that for 7 £ (1/2, 1) the fractionally diffusive (MG) 
equations are locally well-posed in Sobolev spaces H S (T 3 ) for s > 5/2 + (1 — 27). With an additional small 
data assumption we obtain a global well-posedness result. To see why one may not use energy estimates 
to obtain local well-posedness of the (MG) equations when 7 < 1/2, we point out that in the H s energy 
estimate for (1 1 . 10b there are only two terms for which more than "s derivatives" fall on a single function: 

T bad = j A s U-VGA s e and T good = j U • VA S Q A S G (1.12) 

where we denoted A = (—A) 1 / 2 . Since V • U = 0, upon integrating by parts we have T goo d = 0. On 
the other hand, the term T^ad does not vanish in general. The diffusion competing with T\, a d is given by 
«|| A S+7 Q|| 2 2 . Two issues arise. First, since U « AG and 7 < 1/2, and one cannot control || A s ~' y U\\i J 2 
with || A s+7 0|| i 2, and therefore bounding Tb a( i could only be achieved by exploring some extra cancellation. 
Second, since the symbols Mj are even the operator M is not anti-symmetric, thus one cannot re-write Tf, a d 
as — f A S Q • iW[VQ A S Q] (note that if one could do this, a suitable commutator estimate of Coifman-Meyer 
type would close the estimates at the level of Sobolev spaces, as in l]3j[T0l). Hence there is no a commutator 
structure in Tb a( i, and the estimates do not seem to close at the level of Sobolev spaces when < 7 < 1/2. 

In section |4] we consider the case 7 £ (0,1/2) and prove that the solution map associated with the 
Cauchy problem is not Lipschitz continuous with respect to perturbations in the initial data in the topology 
of a certain Sobolev space X. Hence the Cauchy problem is ill-posed in the sense of Hadamard. This is 
achieved by considering a specific steady profile Go, and constructing functions which are "close" to Go, but 
which in arbitrarily short time deviate arbitrarily far from Go, when measured in X. The arguments hinge 
on using techniques of continued fractions to exhibit an unstable eigenvalue for the linearized equation 
(cf. dHHl). These techniques produce a lower bound on the growth rate of eigenvalues to the linearized 
equations, and in the case where 7 £ (0, 1/2) we prove that the magnitude of this lower bound can be made 
arbitrarily large. Once unstable eigenvalues of arbitrarily large part have been obtained for the linearized 
equations, the Lipschitz ill-posednedness of the full non-linear equations may be obtained using classical 
arguments (see, e.g. lfl^[TTl[T2ll20ll24l ). We emphasize that the mechanism producing ill-posedness is not 
merely the order one derivative loss in the map Q — > U. Rather, it is the combination of the derivative loss 
with the anisotropy of the symbol M and the fact that this symbol is even. We note that the even nature 
of the symbol of M plays a central role in the proof of the non-uniqueness given in [22] for the L 00 -weak 
solutions of the non-diffusive (MG) equations. 

In section [5] we study the more subtle transitional case when 7 = 1/2. In this case, if the initial data is 
"small" with respect to k, we use energy estimates to prove that there exists a unique global solution. In 
dramatic contrast, there exists initial data which is "large" with respect to k, and for which the associated 
linear operator has arbitrarily large unstable eigenvalues, which may be use to prove that the equations are 
Hadamard ill-posed. Situations where the qualitative behavior of the solution depends crucially on the norm 
of initial data are also encountered for other evolution equations. For example, in the two-dimensional 
Keller-Segel model of chemotaxis, properties of the solution are strongly dependent on the total mass m 
of cells. If m < 8tt, a global bounded solution of the initial value problem exists and is dispersed for 
t — > 00, while for m > 8tt, blow-up in finite time occurs. The critical case m = 8tt is by now understood 
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and the result states that a global solution exists and possibly becomes unbounded as t — > oo (see lfl3l 
and references therein). We also mention that for the Rayleigh-Taylor and the Helmholtz problems it is 
"geometric" conditions, rather than merely the size of the data, which leads to ill-posedness (see, for instance 
0). 

In section [6] we take advantage of the anisotropy of the symbols dl.4b - dl.6l > and observe an interesting 
phenomenon: when the initial data and the source term are "well-prepared" in an appropriate sense, it is 
possible to prove stronger regularity results for the ensuing solutions. More specifically, when the initial 
data and the force have "thin" Fourier support, i.e. they are supported purely on a plane V q = {k = 
(ki, &2, ks) G : ki = qk\] in frequency space, where q is a fixed non-zero rational number, the operator 
M behaves like an operator of order zero. The reason is that the Fourier symbols Mj are unbounded only 
on curved regions in frequency space (cf. d 1.8b above), but if k G V q the Mj(k) are bounded by a constant 
depending on q. This observation, combined with the fact that when the data and source have frequency 
support on V q , then so does the solution of the (MG) equations at all later times, means that the smoothing 
properties of the fractional diffusion term are stronger than they are for the generic data situation. Hence 
it is possible to prove stronger regularity results, and in particular the local existence and uniqueness of 
Sobolev solutions holds for all values of 7 G (0,1). We note that this is not in contradiction with our results 
proven in sections [4] and [5J since in order to obtain ill-posedness we need to send q — > (so that we obtain 
eigenvalues of the linear operator with arbitrarily large real part), whereas in the case of well-prepared data 
the value of q > is fixed. Other uses of thin Fourier support for different problems can be found, e.g. in 

num. 

2. Preliminaries 

This section contains a few auxiliary results used in the paper. In particular, we recall the, by now 
classical, product and commutator estimates, as well as the Sobolev embedding inequalities. Proofs of these 
results can be found for instance in fl31l23ll24l . Let us denote A = (—A) 1 / 2 . 

Proposition 2.1 (Product estimate). Ifs > 0, then for all /, # G H s n L°°, and we have the estimates 

||A S (/<7)|| LP < C(\\f\\ LPl \\A s g\\ LP2 +||A7|MMIl*4), (2.1) 

where l/p = 1/pi + l/pa = 1/P3 + 1/P4» an d P,P2,P3 £ (lj 00). In particular 

W(fg)\\ L * < C(\\f\\ L ~\\A'g\\v + \\A s f\\ L 4g\\ L oo) . (2.2) 

In the case of a commutator we have the following estimate. 

Proposition 2.2 (Commutator estimate). Suppose that s > and p G (1, 00). If, f,g G S, then 

\\A s (fg) - fA s g\\ LP < C (||V/|| LPl IIA^H^ + ||A7II^3 \\g\\ L ^) (2.3) 

where s > 0, l/p = 1/pi + l/p 2 = 1/P3 + 1/.P4. andp,p 2 ,P3 € (l,oo). 

We shall use as well the following Sobolev inequality. 

Proposition 2.3 (Sobolev inequality). Suppose that q > 1, p G [q, 00) and l/p = 1/q — s/d. Suppose that 
A s f£ L q , then f G L p and there is a constant C > such that 

||/||l p <C||A*/|| L9 . (2.4) 

Throughout this paper we shall make use the following definition of solutions of d 1- lb — (II .3b . 

Definition 2.4 (Solution of the (MG) equation). Let s > 1/2, T > 0, k > 0, and 7 G (0,1). Given 
00 G H S (T 3 ), we call a function 

Q G L°°(0, T; H S (T 3 )) n L 2 (0, T; H s+ ~< (T 3 )) (2.5) 
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a solution of the Cauchy problem (11-1 b — <l 1 - 3 b . if 

f (e(t,-)-Q )(pdx + K f f e(-Ay^dxds 

Jt 3 Jo Jt 3 

holds for all test functions if G C^°(T 3 ), and all t G (0, T). 



JT 3 



?7 • Vtpdxds = 



(2.6) 



3. The regime 1/2 < 7 < 1: well-posedness results 

It has been shown in (9l that the system (ll.lb - dl.6l l is globally well -posed when 7 = 1. In this section 
we prove that in the range 7 G (1/2, 1) equations (II . lb — d 1.6b are locally well-posed in H S (T 3 ). With an 
additional assumption on the size of the initial data and the source term, we obtain a global well-posedness 
result. 

Regarding arbitrarily large initial data, we obtain the following result. 

Theorem 3.1 (Local existence). Let 7 G (1/2, 1), and fix s > 5/2 + (1 - 27). Assume that 6 G H S (T 3 ) 
and S G L°°(0, 00; ii~ s_7 (T 3 )) have zero-mean on T 3 . Then there exists a time T > and a unique smooth 
solution 

9 G L°°(0, T; H S (T 3 )) n L 2 (0, T; i? s+7 (T 3 )) (3.1) 
of the Cauchy problem (11. 10b — dl. lib . 

Proof. We multiply equation dl.lb by A 2s 0, integrate by parts to obtain the following .£P-energy inequality 



^l|A s ef L2 + /,||A^e||i 2 < 



5A 2s e 



+ 



U ■ V6A 2s 6 



We estimate the first term on the right side by 

SA 2s Q 



< \\A S ^S\\ L2 \\A S ^Q\\ L 2 < i-||A s -^||| 2 + |||A s +^e|| 2 2 . 



To handle the second term we proceed as follows. First note that 



U ■ V9A 2s e 



/ a s ^'(u ■ ve)A s+7 e 

JR3 



< \\A s -^u-ve)\\ L 2\\A s ^@\\ L 2. 



(3.2) 



(3.3) 



(3.4) 



The estimate of the product term follows from Proposition ^. II Hence, we have 

||A S ^(?7- ve)|| L 2 < c7(||A s ^c/|| L p||ve|| 2 p +||A s ^ve|| LP ||*y|| 2P ) (3.5) 

LP-' 2 Lp~' 2 

for some 2 < p < oo to be chosen later. Recall that Uj = MjQ, where the symbol of the multiplier Mj 
enjoys a uniform bound |Mf(fe)| < C\k\, for j = 1, 2, 3. Therefore A~ 1 Mj are bounded operators on LP 
(see (9j Section 4] for details) and we obtain 

\\A s ^U\\ L v < C\\A s -^ +1 @\\ LP (3.6) 

and 

(3.7) 



||J7|| 2 P < C\\AQ\\ 2 P 

Lp-' 2 Lp-' 2 

Hence, the product term is bounded by 

\\A S ~^(U ■ V6)|| L 2 < C7||A9n 

We now fix an arbitrary p such that 

3 

< p < 



Lp- 



6 



| A s-7+l0| 



6 



LP- 



(3.8) 



1 



2(1 - 2 7 ) + 3 5 - 4 7 



For our convenience, we choose to work in sub-critical Sobolev spaces. The following proofs can be also earned out in the 
setting of critical Besov spaces _B|.i or in general Bp q . 
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Note that p > 2 since s > 3/2, and the range for p is non-empty since s > 5/2 + (1 — 27). For 7 £ (1/2, 1), 
our choice of p and the Sobolev embedding (Proposition 12.31 ) gives 

|| A ^7+i || LP < ||A s e||^ c ||A s+ TG||^ 2 (3.9) 

where £ £ (0, 1) may be computed explicitly from 7^ = 3/2 — 3/p + 1 — 7. In addition, since Q has zero 
mean, and p > 3/(s — 1), from the Sobolev embedding we obtain 

||AG|| 2 P <C||A s 9|| i2 . (3.10) 

Combining estimates (I3.2b - (I3- 10b gives 

-^||A S G||| 2 + K ||A s+7 e||| 2 < -||A S -^5||| 2 +C||A s G||^||A s+ ^e||^ c (3.11) 

where < C < 1 is as defined earlier. The second term on the right side of ( 13. Ill ) is bounded using the 
e- Young inequality as 

|||A^G||| 2 +C«-W||A-0||^. (3.12) 
and we finally obtain the following estimate 

-^||A S Q||2 2 + ^||A^G||2 2 < -}\k s ^S\\l 2 +CK~^\\k s ®\\^ . (3.13) 

(XT A Ki 

Using Gronwall's inequality, from estimate (13.131 ) we may deduce the existence of a positive time T = 
r(||5|| L co( ,T;H.-7),||©o||ff,K) such that 9 £ L°°(0, T; H S (T 3 )) n L 2 (0,T; H S+ ^(T 3 )). Note that we 
have 2(2 — £)/(l — £) > 2, and hence we may not obtain the global existence of solutions from the energy 
estimate (13- 13b . if the initial data has large H s norm. These a priori estimates can be made formal using a 
standard approximation procedure. We omit further details. □ 

The second main result of this section concerns global well-posedness in case of small initial data. 

Theorem 3.2 (Global existence for small data). Let 7 and S be as in the statement of Theorem IJ.il 
and let Go £ -£P(T 3 ) have zero-mean on T 3 , where s > 5/2 + (1 — 27). There exists a small enough 
constant e > depending on k, such that if ||Go||° 2 ||Go||^ a + ||©o||™2 ll^lli^o oo-H 3 —t) — £ ' wnere 
a =1— (5/2 + 1 — 27) / s, then the unique smooth solution G of the Cauchy problem (11.1 Ob — (1 1 - 1 1 b is global 
in time, i.e. 

G £ L°°(0,oo;H s (T 3 )). (3.14) 
Proof. We proceed as in the proof of Theorem l3.ll The product term in (13.51 ) is now estimated by 

\\A S -^(U- VG)|| L2 < C(||A S -T£/|| LP ||VG|| 2p +\\A s ^Ve\\ LP \\U\\ 2p ), (3.15) 

LP-' 2 Lp-' 2 

where 

6 

P ~ 2(1 -2 7 ) +3 

so that 

\\A S+1 -^Q\\ LP < C||A S+7 G|| L 2. (3.16) 
With this choice of p and the above embedding, the product estimate gives us 

\\A S ^(U- VG)|| L2 < C||VG|| 2P ||A S+ TQ|| L2 . (3.17) 

Combining (13.171 ) with (13.21 ) and proceeding as in (13.131 ) we obtain 

^4||A s G||| 2 +k||A s+ tQ||2 2 < ||A^S|| L2 ||A S+ ^Q|| L2 + C||VG|| 2p \\A S+ ^Q\\ 2 L2 , (3.18) 
2 at lp-' 2 



SUPERCRITICALLY DIFFUSIVE (MG) 



which in turn implies 



J4l|A s e||| 2 + ^||A^e||| 2 < ±\\A S ^S\\ 2 L2 + C||VG|| 2, ||A s+7 e||| 2 . (3.19) 

z at A ZK LP-' 2 



Observe that 

\l-a 



||ve|| i2p / (P - 2) <c||G||2 2 ||A s e||^ a , 0.20) 

where a = I — (5/2 + 1 — 2j)/s. Therefore, if 



||e||2 2 ||A s eni 2 Q <^, (3.2i) 

estimate (13.191 ). combined with the Poincare inequality ||A s+7 0|| i2 > ||A s @|| i 2, shows that 

4l|A*e||! 2 + £||A'e||£ a < 7 . (3.22) 
at Z k u t n * 



and hence 

||A s 9(t, < ||A s e ||i 2 + _ 7 (3.23) 

for all t > 0. Note that taking the L 2 -product of the equation with gives for any t > 

+ K I W K1 ®Wh dr < ||e ||| 2 , (3.24) 
Jo 

which gives us a basic uniform estimate of in Uj°I? x . Hence, from (13.231 ) and (13.241 ) we obtain that 
condition (13.211 ) is satisfied for all t > as long as we have 

||0o||2 2 ||A s e o ||i 2 Q + ||0o|l2 2 ||A s " 7 ^r)llIr (o>oo;i2) <e, (3.25) 
where e is sufficiently small, thereby concluding the proof of the theorem. □ 



4. The regime < 7 < 1/2: ill-posedness results 

We now turn to the situation when 7 G (0, 1/2), k > 0, and prove that there exists an example of an 
initial datum and time independent force for which it is possible to prove that the fractionally diffusive 
(M G) equation is ill-posed in the sense of Hadamard. The arguments by which we prove this follow the 
lines of the proof in the case when k = in iflOl . We sketch here the proof of 7 G (0, 1/2), but as we shall 
prove below in Section 15.11 the result is also true when 7 = 1/2 and the initial data is large in a suitable 
sense. 

4. 1 . Linear ill-posedness in L 2 . The classical approach to Hadamard ill-posedness for non-linear problems 
is to first linearize the equations about a suitable steady state 0o, and then prove that the resulting linear 
operator L has eigenvalues with arbitrarily large real part in the unstable region (see for instance |[T0l [TT1 
[121201, and references therein). Let ©o = asin(mx^) for some positive amplitude a and integer frequency 
m > 1, and define S = mm 2 " 1 sin (771x3). It i s not nar d to verify that O is indeed a steady state of 
(11-1 b — (1 1 - 2 b - We consider the linear evolution of the perturbation 9 = — Qq, obtained from linearizing 
(fTTTT>-(fTT2T> 

d t 9 + am cos (mx 3 )M 3 6* + k(-A) 7 6» = 0, (4.1) 

where we recall that M3 is the Fourier multiplier with symbol defined in (I1.6I ). The following theorem states 
that the linear equation (14.11 ) is Hadamard Lipschitz ill-posed in L 2 . 

Theorem 4.1 (Linear ill-posedness). The Cauchy problem associated to the linear evolution 



d t 9 = L6 



(4.2) 
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where the linear operator L and the steady state @q are given by 

lo(x, t) = -M 3 e(x, t) d 3 e (x 3 ) - K (-Aye (4.3) 

©0(^3) = asm(mx 3 ) (4.4) 

for some a > 0, integer rn > 1, and 7 E (0, 1/2], is ill-posed in the sense of Hadamard over L 2 . More 
precisely, for any T > and any K > 0, there exists a real-analytic initial data 6(0, x) such that the 
Cauchy problem associated to d4.2| )- (|4.4| ) has no solution 9 £ L°°(0, T; L 2 ) satisfying 

sup <iq0(O,O||y (4.5) 

te[o,T) 

where Y is any Sobolev space embedded in L 2 . 

As we mentioned above, in order to prove Theorem 14. II we construct a sequence of eigenvalues a^' of 
the operator L, which diverge to 00 as j — > 00. For this purpose, given any fixed integer j > 1 we seek a 
solution 9 to (14.21 ) of the form 

6(t, x) = e at sin(j 2 xi) sin(jx2) c p sm(mpx 3 ) (4.6) 

p>i 

where the sequence c p decays rapidly as p — > 00. We shall construct such a solution 8 with er real and 
positive, and in addition obtain bounds on the value of a. Inserting (14.61 ) into (14.11 ) and using the definition 
of M 3 , we obtain 

a„ Cp sm(mpx 3 ) + \^ — (sm(m(p + l)x 3 ) + sin(m(p — 1)2:3)) = 0, (4.7) 
P>1 P>1 ' 

where we have denoted 



and 



a p = a + K{j A +j 2 + {mp) 2 Y< (4.8) 
_ 2 3 n 2 (mp) 2 (f + j 2 + (mp) 2 ) + 2p 2 f 

for any integer p > 1 (note that j is fixed). Here fi = (3 2 /rj. An essential feature of the coefficients a p 
is that they grow rapidly with p as p — > 00. Equation ( 14.91 ) gives the recurrence relation for the unknown 
coefficients c p : 

o- p c p + + = 0, for p > 2 (4.10) 

CT1C1 + — = 0, forp = 1. (4.11) 

Note that given any a > (which then uniquely defines <7 p for all p > 1 cf. (14.81 )) and given any c\ > 0, 
one may use the recursion relations (14. 1 0b — (14. Ill) to solve for all c p with p > 2. However, only for suitable 
values of a do the c p 's vanish sufficiently fast so that 9 is C°° smooth. In this direction we prove: 



Lemma 4.2. Let a p be defined by (14.91) . where the positive integers m and j are fixed, and k, p, Q, a are 
fixed physical parameters. Also, define c\ = a±. Then there exists a real eigenvalue a = a^' > 0, and 
rapidly decaying sequence {c p } p >2 which satisfies (14- 10b — (14- lib - Furthermore the lower bound 

a U) > a^mj 2 (f+j 2 ) ,-4 + , 2 + 4m 2yr (4 12) 

> 2^ 2 m 2 {f + j 2 + 4m 2 ) + 2p 2 f U + 3 + > K ] 

holds, andc p = 0(C p /(p— l)! 4 ) asp — > 00, for some constant C which depends on j and the physical 
parameters. 
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Proof of Lemma WT[ We define 



and write (14. 10b — d47TTT ) in the form 



ri P = C -^^ (4.13) 

Cp—l Olp 

1 

o-pcep + rj p+1 + — = 0, p > 2 (4.14) 
Vp 

(Tlttl + 7/2 = 0, p=l. (4.15) 

Using (14.141 ) to write 7/2 as an infinite continued fraction and equating with 7/2 given by (14.151 ) gives the 
characteristic equation 

<j\ol\ = 3 . (4.16) 

£j 2 a2 1 

T3«3 i 

Recalling the definition of a p for p > 1 cf. (I4.8I ). we observe that (14.161 ) is an equation with only one 
unknown, namely a. For real values of a we define the infinite continued fraction F p (a) by 

F p {a) = l — l (4.17) 



P U P „ 1 

^+l^+l- CTp+2>+2 _ 



and the function 



a p a p - Ja^aj - 4 



G P {<?) = \ = , (4.18) 

a p a p + Ja^aj - 4 

for all p > 2, and all cr such that 02OL2 > 2. We note that due to the very rapid growth of the a p s, for real a 
the continued fraction F p defines a function which is smooth except for a set of points such with 0202 < 2. 
For the rest of the proof we will always assume that a is real such that cr2«2 > 2, which also implies that 
a p a p > 2 for all p > 2. 

We note that by construction G p satisfies 

G P (a) = = —I • (419) 

a p a p - G p {a ) a v a p - 1 — 

P P CTpCtp — ... 

Since we have a p a p — > 00 as p — > 00, we have that G p (a) — > as p — > 00 for every fixed a. Also, we 
clearly have 

G2O) > G 3 (a) > G 4 (a) > . . . > (4.20) 

and 

G p+1 (a) < G p (a) < a p a p (4.21) 
for all p > 2. Hence, from ( 14.191 ) and (14.201 ) we obtain 

G 2 (a) > l — — > l — > 0. (4.22) 

a 2 a 2 - G 3 (a) a 2 a 2 - azaz _ Gi{a) 

An inductive argument then gives 

G 2 {a) > — j = F 2 (a) > (4.23) 

o 2 a 2 1 

J (7404 — ... 

for all a such that 17202 > 2. Repeating this constructive argument we obtain that ^3(0") > ^3(0") > 0, and 
hence 

< a 2 a 2 - Gz{a) < a 2 a 2 - F 3 (a) < a 2 a 2 (4.24) 
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so that 

F 2 {a) = l —— > — . (4.25) 

G2OL2 — r 3 {a) a<i<J2 

In summary, we have proven that 1/(17202) < F<i(p) < ^2(0") for all a such that <72a2 > 2, and the 
straight line aicti intersects both the graph of 1/ ((T2&2) an(1 the graph of (^(cr) in this rage of a. Hence, it 
follows that there exists a real solution a^' of the equation 

F 2 (a (J) ) = Jpa u (4-26) 

with a^a2 > 2, where we denote a p = + k(j 4 + j 2 + m 2 p 2 ) 7 , for all p > 2. That is, is a real 
positive solution of the characteristic equation (14.161) . Furthermore, due to (14.231) and (14.251) after a short 
calculation we obtain an upper and a lower bound on namely 

< a 0) a U) < (4 .27) 



CilOt2 OL\OL2 

For <tw) satisfying ( 14.26b , we now construct the sequence c p which decays rapidly as p — > 00. The 
recursion relation ( |4.14| )- (|4.15l ) uniquely defines the values of rj p . After letting c\ = a\, we define 

c p = a p ??p??p_i ...r]2 (4.28) 

for all p > 2. This sequence satisfies d4- 10b — (14- lib by construction. Furthermore, we observe that rj p = 
—F p (a^>). Repeating the arguments which gave ( 14.231 ) and ( 14.251 ) with 2 replaced by p gives 

< Vp < (4.29) 



«(j 4 + j 2 + 4m 2 y < a (j) (4.30) 



U) lp U) 

(Tp Oip C7 p Up 

Moreover, from ( 14.91 ) we have that a p = 0{p A ) as p — > 00, and hence from ( 14.291 ) we obtain that rj p = 
0(p~ 4 ). Thus it follows from (l4~28l that c p ->■ as p — > 00, and this convergence is very fast, namely 
0(C p (p — l)! -4 ) as p — > 00, for some positive constant C = C(o~v\/j,, a, m,j). Therefore the solution 
given by ( 14.61 ) with a replaced by and c p as defined by ( 14.281 ), lies in any Sobolev space, it is 
C°° smooth, and even real-analytic. 

We substitute for a\ and «2 from (14.91 ) into the bound (14.27b to obtain estimates on given by 

afJ-mj 2 (j 4 + j 2 ) 

2 3 ft 2 m 2 (j 4 + j 2 + 4m 2 ) + 2p 2 j 4 
and 

(j) 2a^mj 2 {f + j 2 ) 4 . 2 2 

CT < 0302 2f -4 , -2 , 2^.,2'4 " K C? + .7 + 77J ) . (4.31) 

We recall the role of the constants in the above expressions (14.30b — (14.3 1 b for the bounds on The 
physical parameters are k the coefficient of thermal diffusivity, Q the rotation rate of the system, /j, is related 
to the underlying magnetic field, and a is the magnitude of the steady buoyancy @o- Expressions (14.30b 
and (14.31b indicate wide ranges of the physical parameters for which there exist eigenvalues which are 
positive as postulated in the construction of the solution to (14.1b . We note that in the Earth's fluid core k is 
very small. □ 

Lemma 4.3. Let 7 £ (0,1/2). Fix a,p,,Q,K > and an integer m > 1. There exists an integer 
jo = jo( a i m i ^> K ) 7) such that for all j > jo there exists a C°° smooth initial datum 6^(0, x) with 
\\6^ (0, -)||l 2 = 1 an d a C°° smooth solution 9^\t,x) of the initial value problem associated with the 
linearized (MG) equation (14.2b . such that 

\\e^\t,-)\\ L 2 >exp(j 2 CU) (4.32) 

for all t > 0, where C* = C* (a, m, fj,, K, jo) is a positive constant defined in (14.33b below. 
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Proof of Lemma W3\ For any j > 1, Lemma l431 guarantees the existence of an eigenvalue irW of the linear 
(MG) operator L, with associated C°° smooth eigenfunction 9^(0,x) given by letting t = in (14.6I ). 
Then6>^(t,x) = exp(<7^t)6^(0, •) is a solution of (I4.2I ). It is clear that one may normalize 9^'(0, ■) to 
have L 2 norm equal to 1, and hence ||#W(t, -)\\ L 2 = exp(cr^H). 

The lemma is then proven if we pick j large enough so that a^>> > j 2 C* for some positive constant C* 
(independent of j). This is guaranteed by the lower bound (14.301 ) for if we pick j > jo, where jo is a 
large enough fixed integer such that 

Note that when 76 (0, 1/2) such a jo always exists, independently of the size of the physical parameters. 

□ 

We now have all necessary ingredients to conclude the proof of Theorem 14. II 

Proof of Theorem WA\ Let T > and K > be arbitrary, and let Y C L 2 be a Sobolev space. Pick an 
integer m > 1, and an amplitude a > 0. For these fixed a, m and physical parameters fi, Q, k, 7 > 0, define 
jo and C* as in Lemma 1431 For any j > jo, Lemma 1431 guarantees that there exists a C°° smooth initial 
condition 6^>\0,x) which we re-normalize to have || (0, -)||y = 1, such that the associated solution 
9^\t, x) of the Cauchy problem (H~2l-(14~I41 satisfies 

\\e^{t,-)\\ L , > exp(j 2 C*t)\\0 {j) (0,-)\\ L2 (4.34) 

for all t > 0. We note that this solution to the linear equation is the uniqud3 solution in L°°(0, T; L 2 ). We 
now claim that there exists a sufficiently large > jo sucn that 

\\e^\T/2,-)\\ L 2>2K (4.35) 

which would then conclude the proof of the Theorem, since (0, -)||y = 1. After a short calculation it 
is clear that (14.351 ) follows from (14.341) if we manage to prove that 



cxp 



') ||0 W (O, >2K. (4.36) 



But from the definition of 0(>) cf. $M above, we see that ||6K>)(0, -)|| L2 > c x = a\ > l/(Cj 2 ), for some 
sufficiently large constant C, cf. ( 14.91 ) above. The fact that for every fixed T > and C* > the sequence 
exp(j 2 C*T/2)/(Cj 2 ) diverges as j — > 00, shows that there exits some sufficiently large j* such that (14.361 ) 
holds, thereby concluding the proof of the theorem. □ 

4.2. Non-linear ill-posedness in Sobolev spaces. Having established that the linearized (MG) are ill- 
posed in L 2 , we now turn to address the Hadamard ill-posedness of the full nonlinear (MG) equations. The 
ill-posedness of non-linear partial differential equations may have different sources, of varying degree of 
gravity: finite time blow-up, non-uniqueness, or non-smoothness of the solution map, to name a few. As in 
lllOl for the case k = 0, here we show that the fractionally diffusive (MG) equations, with < 7 < 1/2 
do not posses a Lipschitz continuous solution map. Recall the definition of Lipschitz well-posedness for the 
nonlinear equation (see lfT2l Definition 1.1], or lITOl Definition 2.9]): 

Definition 4.4 (Lipschitz local well-posedness). Let Y C X C W 1,4 be Banach spaces. The Cauchy 
problem for the (MG) equation 

d t e + u • ve + K(-A) 7 e = s (4.37) 

V • U = 0, Uj = MjO (4.38) 



n 

' The proof of uniqueness of finite energy solutions to the linear equation follows from a representation of the solution as a 
Fourier series (see 1101 . Proposition 2.8). 
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is locally Lipschitz. (X, Y) well-posed, if there exist continuous functions T, K : [0, oo) 2 — > (0, oo), the time 
of existence and the Lipschitz constant, so that for every pair of initial data 6^ (0, •), (0, •) G Y there 
exist unique solutions Q^,Q^ G L°°(0,T;X) of the initial value problem associated to (14.37 b — (14.38 b . 
that additionally satisfy 

\\@ {1) (t, •) - e {2 \t, oiu < ^l|e (1) (o, •) - e(2) (o, -)\\ Y (4.39) 

for every t e [0,2], wfereT = T(ll© (1) (0, Oily, ll@ (2) (0, -)\\y) and K = K(\\@W(0, -)\\ Y , \0 2 \o, -)||y). 

Remark 4.5. Clearly the time of existence T and the Lipschitz constant K also depend on ||5||x,oo(o,co;y)> 
and on n, but we have omitted this dependence in Definition 14 .4 1 since it is the same for both solutions 
and 6^. 

Remark 4.6. If e^(t,-) = and X = Y, Definition S3 recovers the usual definition of local well- 
posedness with a continuous solution map. However, Defintion l4.4l allows the solution map to lose regularity, 
which is usually needed in order to obtain Lipschitz continuity of the solution map for first order quasi-linear 
equations. Hence, the typical pairs of spaces (X, Y) that we have in mind here are X = H s , and Y = H s+1 , 
with s > 1 + 3/4. 

The main result of this section is the following theorem. 

Theorem 4.7 (Nonlinear ill-posedness in Sobolev spaces). The (MG) equations (14.37b - (14.38b . with 7 G 
(0, 1/2) are locally Lipschitz (X,Y) ill-posed in Sobolev spaces Y C X embedded in W 1A (T 3 ), in the 
sense of Definition \4.4\ 

For the purpose of our ill-posedness result, we shall let ®( 2 \t,x) be the steady state 6(0:3) intro- 
duced earlier in (14.4b . We consider X to be a Sobolev space with high enough regularity so that dt@ G 
L°°(0,T; L 2 ), which implies that G is weakly continuous on [0, T] with values in X, making sense of the 
initial value problem associated to (14.37b — (14.38b . The proof of Theorem 14.71 follows from the strong linear 
ill-posedness obtained in Theorem [4j] and a fairly generic perturbative argument (cf. |[20l Thorem 2] or |[24l 
pp. 183]). The proof of Theorem 14.71 follows the lines of the proof for the non diffusive problem given in 
HQ). 

Proof of Theorem W7\ Since the Sobolev space X embeds in H l , and 7 G (0, 1/2), the linearized operator 
LQ = — M%® d%®Q — k(— A) 7 © maps X continuously into L 2 , and since X C W 1,A , the nonlinearity 
Are = -MjOdj® is bounded as ||JVe|| £ 2 < ||V9|| 2 4 < C||0|&, for some constant C > 0. 

Fix the steady state ©0(^3) G Y, as given by (14.4b . Also, fix a smooth function ^0 £ Y, normal- 
ized to have ||^o||y = 1, to be chosen precisely later. The proof is by contradiction. Assume that the 
Cauchy problem for the (MG) equation (14.37b — (14.38b is Lipschitz locally well-posed in (X, Y). Consider 
9( 2 )(0,a:) = 60(3:3), so that Q^ 2 \t,x) = 6 (x 3 ) for any t > 0. Also let 

6 e (0,a:) = eo(x3) + eVo(zO, 

for every < e < || ©0 II V- To simplify notation we write 6 e instead of 6^ 1,e ^. By Definition 14.41 for 
every e as before there exists a positive time T = T(||6o||y, ||6 e ||y) and a positive Lipschitz constant 
K = (||6o||y, ||6 e ||y) such that by (14.39b and the choice of ifio we have 

||6 e (t,-)-6o(-)||x <Ke (4.40) 

for all t G [0, T\. We note that since ||6 e (0, -)\\y < \\&o\\y + e < 2||6 ||y, due to the continuity of T and 
K with respect to the second coordinate, we may choose K = i^(||6o||y) > and T = T(||6o||y) > 
independent of e G (0, ||6 ||y), such that (14.40b holds on [0, T]. 
Writing 6 e as an 0(e) perturbation of 60, we define 

nt,x) = @e(t,X) ~ @0{X3 \ (4-41) 
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for all t € [0, T] and all e as before. It follows from (14.401 ) that ip e is uniformly bounded with respect to e 

in L°°(0, T; X). Therefore, there exists a function ip, the weak-* limit of ip € in L°°(0, T; X). Note that V> e 
solves the Cauchy problem 

d t tp e = Lip e + eN(ip e ) (4.42) 

^(0,-)=^ - (4-43) 

Due to the choice of X, we have the bound 

\\Nr\\Li<C\\^\\ 2 x <CK 2 , (4.44) 

and from (14.421 ) we obtain that dtip e is uniformly bounded with respect to e in L°° (0, T;L 2 ). Therefore the 
convergence ip e —> ip is strong when measured in L 2 . Sending e to in (14.421 ). and using (14.441 ). it follows 
that 

d t ip = Lip 
^(0, •) = ipo 

holds in L°° (0, T; L 2 ), and this solution is unique since the problem is now linear. In addition, the solution 
ip inherits from (14.401 ) the upper bound 

\W,-)\\v<K (4.45) 

for all t € [0, T]. But this is a contradiction with Theorem 14. II Due to the existence of eigenfunctions for 
the linearized operator with arbitrarily large eigenvalues, one may choose ipo (as in Lemma 1431 ) to yield a 
large enough eigenvalue so that in time T/2 the solution grows to have L 2 norm larger than 2K (similarly 
to the proof of Theorem 14 .II ). therefore contradicting (14.45b - □ 

5. The regime 7 = 1/2: a dichotomy in terms of the size of the data 

If the initial data is small with respect to k, then one may use energy estimates to show that there exists a 
unique global smooth solution evolving from this data (cf. Section l5T2"l below). However, the proof does not 
apply for the case of large data, not even to prove the local existence of solutions. In the case of large data, 
we may construct a steady state such that solutions evolving from initial data which is arbitrarily close to 
this steady state diverge from it at an arbitrarily large rate, for positive time (cf. Section IBTTl below), i.e. the 
equations are Hadamard ill-posed. 

5.1. IU-posedness for 7 = 1/2 and large data. In Section [4] above we have shown that for 7 < 1/2 the 

(MG) equations are Hadamard ill-posed in Sobolev spaces, in the sense that there is no Lipschitz continuous 
solution map (see Theorem I4.7I ). The main ingredient in the proof of ill-posedness for 7 £ (0, 1/2) was to 
show that the linearized (MG) operator 

L9 = -M 3 6d 3 G - K(-Aye (5.1) 

where 

0o = asin(mx3) (5.2) 

is a steady state associated to the source term S = nam 2 " 1 sin(mx3), has a sequence of eigenvalues a^' 
which diverge to infinity as j — > 00 (see Lemmas l4.2l and 1431 ) . 

We emphasize that it is only in the proof of Lemma 1431 where 7 < 1/2, rather than 7 < 1/2 is used. 
Indeed, in order to prove that given any real a > and any integer m > 1 there exists some sufficiently large 
integer jo such that the constant C* of (14.331 ) is strictly positive, 7 < 1/2 is both necessary and sufficient. 
In the case 7 = 1/2, we can prove the following alternative to Lemma 1431 which states that only if a is 
sufficiently large with respect to k (in terms of \i and O), then one obtains a sequence of eigenvalues which 
diverge to 00. Namely: 
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Lemma 5.1. Let 7 = 1/2. Fix a, fi, 0, and m. For values of a such that 

K < 25ft2 m 2 + 2 /i 2 ( 3) 

and all integers j > m, the statement of Lemma \4. 3\ holds with the constant C* given by 

_ ap,m 

°* ~ 25Q2 m 2 + 2fJ 2 ~ K - V> A > 

Proof of Lemma I3T71 The proof follows from the lower bound ( 14.301 ) on the eigenvalue similarly to the 
proof of Lemma [4731 The role of condition ( 15.31 ) is to guarantee that there exists a large enough j such that 

q/xm(j 4 +j 2 ) «_ 4 2 2 

2302 m 2 (i 4 + Jo 2 + 4m 2) + 2/i2j 4 ^ Uo +^0+^) > U. 

To avoid repetition we omit further details. □ 

For those values of a for which (15.31 ) holds, Lemma lBTTl shows that the operator L has unbounded spectrum 
in the unstable region, and hence one may prove with virtually no modifications from the 7 G (0, 1/2) case, 
that the equations are ill-posed. More precisely we have 

Theorem 5.2 (Ill-posedness for large data). LetY cX c W 1,4 (T 3 ) be Sobolev spaces, and let 7 = 1/2. 
Given k, > 0, fix an integer m > 1. Le? a > Z?e sufficiently large such that (15.31 ) holds and let 



Oq 1 ^ = asin(mx3) £ Y be a steady state of d 1 - lb - Zfoen, g/verc a^ry T > a7i<i anj if > 0, f/iere ex/sfs 
an initial condition G y amc? a corresponding solution 0( 2 ) G L°°(0,T;X) o/f/je Cauchy problem 

(fTTTb— (fTT3T>. swc/i ^a/ 

sup ||e> (2) (i, -) -e«(t,-)IU > 2K||eJ 2) -e«||y (5.5) 

te[o,r] 

where Q^(t, •) = Gq 1 ^. 77iaf jj, f/ze (MG) equations are locally Lipschitz (X,Y) ill-posed when the data is 
large. 

The proof of Theorem l5.2l is the same as the proof of Theorem 14.7 1 above. The only difference in the case 
7 = 1/2 is to use Lemma lBTTI to show that the linear equations have unbounded unstable spectrum. We omit 
details. 

5.2. Well-posedness for 7 = 1/2 and small data. As shown in the previous section, we can exhibit initial 
data, for which the system (II .IN 1.21 ) is ill-posed in the above described sense. Note also, that the proof of 
Theorem 13. II fails for the value 7 = 1/2. Thus, 7 = 1/2 indeed is the limit of the local well-posedness 
theory. Nonetheless, we still can prove that the considered system is globally well-posed for small data. 

Theorem 5.3 (Global existence for small data). Lets > 5/2 and assume that the initial data Qq G H s (T 3 ) 
and S G L 2 (0, 00; -fP _7 (T 3 )) have zero-mean on T 3 . There exists a sufficiently small constant e > 
depending on k, such that if WGoWlzWQoW 1 ^ + ||5|| L°°(o,oo;_ff a -T) — e > where a = 1 — (5/2)/ s, then the 
unique smooth solution 

6 G L°°(0,T;H s (T 3 )) (5.6) 
of the Cauchy problem (11.11 )— d 1.61 ) is global in time. 

Proof. We proceed as in the proof of Theorem [3^2] and obtain the energy estimate 



^||A S 6||| 2 + k||A s+ ^9|| 2 2 < ||A s -ts|| L 2||A s+ ie|| L 2 + 



/ k s ~^(U- VG)A s+ 29 



(5.7) 



The second term on the right side is estimated using the product estimate in Proposition 12. 11 Thus we obtain, 
similarly to ( 13.81 ) 

\\k s ^(u ■ ve)|| L2 < c(||ve|| L cc + ||c/|| LOO )||A s +5e|| L2 . (5.8) 
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Since s > 5/2, we get 

\\u\\ L oo + ||ve|| L ~ < c||e||22||A s e||^ Q , (5.9) 

where a = 1 — (n/2 + l)/s. Combining estimates (I5.7I )- (|5.9I ) gives 

l^M s ®\\h + < ^\\^s\\h + ^||A s +ie||| 2 + ciieii^iiA-eifciA'+^eiii,, 

Z at Zk Z 

(5.10) 

which in turn leads to 

i^||^V-0||i 2 < ^NyV—^S'lli, ^ CIIOII^^II^V'OII^^II^V'^^QIlis- (5-11) 

We obtain the desired result as in the proof of Theorem 13.21 □ 

Remark 5.4. Note that conditions in Theorem |3.2| are consistent with the above theorem if we set 7 = 1/2. 

Remark 5.5. We note that the ill-posedness result in the case 7 = 1/2 requires the constant defined in 
(15.41) to be positive. This does not hold when the value of a is sufficiently small (depending on k, //, and 
H) no matter what value of m > 1 is picked. Recalling that a measures the magnitude of the initial data 
associated with ||Oo||z,2> we observe that the well-posedness result when 7 = 1/2 is only proven in the case 
of small data and small force (see Theorem 15 .3 1 ), and hence for a sufficiently small. Therefore our large data 
ill-posedness result is consistent with the small-data well-posedness result when 7 = 1/2. 

6. Improvement in regularity for "well-prepared" data and source 

In this section we explore the following observation: if the frequency support of G lies on a suitable 
plane, then the operator M is "mild" when it acts on Q, i.e. it behaves like an order operator, and hence 
the corresponding velocity U is as smooth as G. This enables us to obtain improved well-posedness results 
over the generic setting when no conditions on the Fourier spectrum of the initial data (and source term) are 
imposed. For instance the local existence an uniqueness of smooth solutions holds even if < 7 < 1/2, a 
setting in which we know that for generic initial data the problem is ill-posed. 

To be more precise let q G Q be a rational numbei0 with q 7^ 0. We define the frequency plane V q to be 
the set 

V q = {k = (h,k 2 ,k 3 ) £Z 3 :k 2 = qh}. (6.1) 
We shall need the following straightforward observation: 

Proposition 6.1. Assume f, g are smooth T 3 periodic functions, with frequency support supp(/), supp(g) C 
Vqfor some q € Q \ {0}. Then we have supp(/ ± g), supp(/g), supp(A-lj/) C V q ,for all j € {1, 2, 3}. 

Proof of Propostion \6. 1 1 Clearly supp(/ + g) = supp(/ + g) C supp(/) U supp(<?) C V q . To see that 
that the frequency support of the function Mjf lies on V q , we just note that supp(Mjg) = supp(Mj/) C 
supp(Mj) n supp(/) C V q . To analyze the frequency support of fg = f * g, note that supp(/ * g) C 
supp(/) + supp(^) C V q , since V q is closed under addition of vectors, concluding the proof. □ 

Lemma 6.2. If f is a smooth T 3 periodic function with f(k\, k2,0) = Ofor all k\,k 2 £ Z, and such that 
supp(/) C Vqfor some q € Q\ {0}, then there exists a universal constant C* = C*(q,£l, /3,r]) > such 
that 

(Mj~f)(k)\ = \Mj(k)f(k)\ < C*|/(fc)| (6-2) 
for all k £ Z 3 , and for all j G {1, 2, 3}. Additionally, the constant C* blows-up as q — > 0. 



Since we are in the periodic setting when the frequency is supported on Z . 
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Proof of Lemma \6?2\ It is clear that (16.21 ) has to be proven only for k such that ks 7^ 0, since otherwise we 
have that f(k) = and the statement holds trivially. We now consider each of the cases j G {1, 2, 3}. 
Without loss of generality take q > 0. 

For j = 1, a short algebraic computations gives 



Mi(fc; 



from which it follows that |Mi(fc)| < C, for a suitable constant C, by using the inequality k\k% + fci&l < 
fe| + fc|. Similarly to (16.31) it follows that for fe 6 P 9 we have 

awl^w^^w ^w ,6 ' 4) 

and 

M 3 (*)| < (, 2 d + rt^) <6 - 5> 

which concludes the proof of the lemma upon using kfk^ + k\k\ < k\ + fc|. □ 

6.1. Local existence and uniqueness for < 7 < 1. The main result of this section is: 

Theorem 6.3 (Local existence). Let 7 e (0, 1), and fix s > 5/2 — 7. Assume ?/jct? Go € -£P(T 3 ) and 
S € L°°(0, 00; .fP _7 (T 3 )) have zero-mean on T 3 anc? are such that 

supp(6o) C V q and supp(5(t, •)) C V q , (6.6) 

for some fixed q £ <Q> \ {0} a«J all t > 0. 77ie72 ?/zere exisft? a T > <3«<f a unique smooth solution 

6 £ L°°(0, T; H S (T 3 )) n L 2 (0, T; ir +7 (T 3 )) (6.7) 

of the Cauchy problem (11.1 0b — d 1.11b . such that 

supp(G(t, •)) C V q (6.8) 

for all t e [0, T). 

Proof of Theorem lO] In order to construct the local in time solution ©, with frequency support on V q , 
consider the sequence of approximations {0 n } n >i given by the solutions of 

tQl + (-A) 7 6i = S (6.9) 

6i(0,-) = e (6.10) 

and respectively 

d t Q n + c/ n _i • ve n + (-A) 7 e n = s (6.11) 

(7 n _i = M6 n _! (6.12) 

e„(o,-) = e (6.13) 

for all n > 2. One may solve (16.9b - (l6.10b explicitly using Duhamel's formula, and noting that the operator 
exp((— A) 7 i) is a Fourier multiplier with non-zero symbol given by exp(|fc| 7 i), it follows from Proposi- 
tion l6.1l that supp(Oi(i, •)) C V q for all t > 0. In addition, we have the bound 

||A S Gi|| Ltx .(o,T;Z, 2 ) + II^ S+7 ©iIIl2(0,T;L 2 ) - II^©o|Il2 + r||5|| ioo ( 0)T ;H s —<) ( 6 -14) 

for all T > 0. On the other hand, in order to solve (I6.11b - (l6.13b we appeal to Theorem IA.1I Indeed, 
by the inductive assumption we have that & n -\ € L°°(0,T; H s ) and also supp(O n _i(T, •)) C V q for 
all T > 0. Hence, by applying Lemma I6T21 we have U n -\ € L°°(0,T; H s ) and by Proposition 16. 1 1 we 
have supp(L 7 " n _i(i, •)) for all t > 0. Therefore all the conditions of Theorem |A. 1 1 are satisfied, by letting 
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v = U n -!, and there hence exists a unique solution 9 n G L°°(0, T; 7P)nL 2 (0, T; H s+ ^) of (l67TTb - (f67T3l> . 
such that supp(9^(t, •)) C V q for all t € [0, T). 

To prove that the sequence {© n } converges, we first prove that it is bounded. Fix a time T to be chosen 
later such that T < || A S O ||^ 2 /||5'||| OO(O T . H s--,y It hence follows from (16. 141 ) that 

IIA^Hl^o,^) + HA^ejl^^ < 2||A s e ||! 2 (6.15) 

holds when j = 1. We assume inductively that (16.151 ) holds for all 1 < j < n — 1, and proceed to prove that 
it holds for j = n. Since U n -i is obtained from Q n -i by a bounded Fourier multiplier (cf. Lemma |6T2~1 ). 
there exists a positive constant C q , depending on q (and on the physical parameters Q,, /3, rf), such that 

\\A s U n ^(t, OHia < C g ||A'e^i(t, Ollia (6-16) 

for all f > 0. Applying A s to (16.111 ) and taking an L 2 inner product with A s Q n we hence obtain 

^IIA^eji. + iiA^iii. 

< ||A a e n || L3 1| [I/^i • V,A s ]@ n \\ L2 + ||A»-^e n || £2 ||A--T5'|| i a 

< c||A s e n || L2 (||vc/ n _i|| L 3/ 7 ||A s e re || i6 /(6- 7 ) + HA'i/n-iii^iiveniiLco) + ||A^e n ,|| L 2||A s -^|| i2 

< C||A s e n || i2 ||A s C/ n _ 1 || i2 ||A s+ ^e ri || L2 + ||A s+ T0 n || i2 ||A s -^|| i2 

< l||A s +^e„||| 2 + C' s , 7 ||A s J7 n _ 1 ||| 2 ||A s G n ||| 2 + ||A s ^5||| 2 (6.17) 

where C sn > is a constant. Above we have used the fact that V • C7 n _i = in order to write the 
commutator, and also the product and commutator estimates of Proposition 12.11 Using the bound (16. 161 ). 
and the inductive assumption (I6.15I ). it follows from (16.171 ) that 

^||A s 9 n ||2 2 + ||A s+ ^e n ||i 2 < 4C S)7 a g ||A s eo||| 2 ||A s G„||i 2 + 2||A S -^||| 2 (6.18) 

and hence 

l|A*e n (v)lli 2 + /V^e-Mlli^ 

Jo 

< ||A s e o ||| 2 e 4a - C « l,AS0o|l ^* + re^-^II^II^^IIA'-TSCr, Oilier 

Jo 

< e 4C ^ Ase <^ (||A s e ||| 2 + tWSWl-^Hs-^) (6-19) 
for all t € (0,T). Therefore, letting 

T< . f l|A s 9o||| 2 In 4/3 ) 

1 < mm < — - — -o , -TT- > (6.20) 

" 1 211^11^(0^-.) '4C Si7 C g ||A^G ||i 2 / 

we obtain that from (16.191 ) that (16.151 ) holds for j = n, and so by induction it holds for all j > 1. This shows 
that the sequence 6„ is uniformly bounded in L°°(0, T; H s ) n L 2 (0, T; iF+f). 

Moreover, we may show that the sequence {9 n } is Cauchy in L°°(0, T; 77 s_1+7 ). To see this, consider 
the difference of two iterates Q n := Q n — Q n -i- It follows from (16.1 II ) that @ n is a solution of 

d t e n + (-A)Te n + J7„_i • ve n + • ve„_i = o (6.21) 

9 n (0, •) = (6.22) 
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for all n > 2, where we also denoted U n = MQ n . Applying A s ~ 1+7 to (16.211) . taking an L 2 inner product 
with A S ~ 1+7 Q„, using V • U n ~\ = 0, and the calculus inequalities of Proposition 12.11 we obtain 

^l|A s " 1+7 e n ,lH 2 + ||A-^e n ||i 2 

< ||A s - 1+7 e n || i2 1| [E/ n _! • v,A*- 1+7 ]e n || L2 + iiA-^e^ii^HA- 1 -^^-!. ve n )|| i2 

< c s ||A s - 1+7 e n || L2 (iivc/n-ilUociiA 8 - 1 ^!!^ + HA'-^E/^iHieiivenlUa" 

+ c s ||A s - 1+ Te n || L2 (||ir^ 1 || i =o||A^e ft ||^ + ||A»- 1 -^i7 n _ 1 || L3 ||ve„||L<. 

< C s , g ||A s - 1+7 9 n || L2 (llA^Gn-ill^HA 8 - 1 ^!!^ + ||A s - 1+7 e n _i|| i2 ||A s+7 G n || L2 ) . 

(6.23) 

In the last inequality above we have used the assumption s > 5/2 — 7, and the estimate (16.16I ). It hence 
follows from (16.231 ). (I6.15I ). and the Cauchy-Schwartz inequality that 

sup ||A*- 1+7 e n (t,-)||L 2 < sup \\A s - 1+ ^e n ^(t,-)\\ L 2 /' T ||A s+7 G ri (i,-)||L^U T| i Ai,+ ' T0 "- 1 ( r '- )ll ^ dT )^ 

te[o,T] te[o,T] Jo 

< sup ||A s - 1+7 Q n _ 1 (t,O||L 2 ^C s , g ||A s Qo|| L2 e^ c - l|AS0oll ^dt. 
te[o,T] 



,91 

(6.24) 



Recalling (16.201 ). if we let T be such that 



n^Qollj. ln4/3 

1 = mm < — - — -7: , -7T- , ■ > (6.25) 

2 ll 5 lli- (0iOO;HS - 7) iC s ^C q \\A s @ Q \\ 2 L 2 8exp(^2hU/3C Stq /^/C^ q 



|2 ' An n II As<a_l|2 ' 

we obtain from (16.231) that 



sup ||A s - 1+7 G n ,(V)||L 2 < l sup ||A s_1+7 Q re _i(t, -)lk 2 - (6-26) 
*e[o,T] 1 t&[0,T] 

Thus 6 n is Cauchy in L°°(0, T; # s ~ 1+7 ), and hence G„ converges strongly to G in L°°(0, T; i? s ~ 1+7 ). 
Noting that s — 1 + 7 > 3/2, this shows that the strong convergence occurs in a Holder space rela- 
tive to x, which is sufficient to prove that the limiting function G G L°°(0, T; H s ) n L 2 (0, T; i? s+7 ) n 
Lip(0, T; flTOin(*-27,«-i)) is a solution of the initial value problem JO-JO. 

To conclude the proof of the theorem, we note that if G^ and G^ 2 ) are two solutions of dl.llHll.21 ). then 
G = GW - G( 2 ) solves 

d t G + (-A) 7 Q + C/ (1) • VG + U • VQ (2) = (6.27) 

with initial condition Q(0, •) = 0. An I? estimate on d6.27| ) shows that Q(t, •) = for all t G [0,i), since 
9(2) G L°°(0,T;H s ) and since if G^ 1 - 1 and Q^ 2 ) have frequency support on V q , so does G. □ 

Corollary 6.4 (Local existence for well-prepared data). Let ji,j2 G Z \ {0}, and fix s, 7 as in the 

statement of Theorem \6.3\ Assume that Go G H S (T 3 ) is such that 60(^1,^2^3) 7^ if and only if 
(ki, k-z) = ±(ii, 32) ■ Similarly, assume that forcing S G L°°(0, T s ; H S ~' Y ) is such that S(ki, k2, k%,t) ^ 
if and only if (k± , ^2) = =t(ji j J2), /or all t G [0, T s ). 77ie« f/iere T G (0, T s ] and a unique solution 

G G L°° (0, T; if s ) n £ 2 (0, T; H s+1 ) of the initial value problem dl.llHll.2l ), anJ in addition we have that 
G(/ci, fe, &3,t) = whenever fe/fci / j'2/ii. 

Proof of Corollary \6.4\ Let g = j'2/ii £ Q- The conditions on the initial data and the force imply in 
particular that their frequency support likes on the plane V q . The existence of a unique smooth solution G, 
with frequency support lying on V q follows directly from Theorem l6.3l But {k\, A^, ^3) G V q is equivalent 

to k 2 = qh = 32k\/jx, so that @(ki,k 2 , k 3 ) / only if k 2 /h = j 2 /ji- □ 
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6.2. Global existence and uniqueness for well-prepared data. In the previous section we have proven 
that for 7 G (0, 1), given initial data ©o and source term S which are well-prepared, i.e. they have Fourier 
support on a plane V q for some q > 0, then there exists a local in time solution G of the Cauchy problem 
(11-1 b — (1 1 - 3 b - which has the property that its Fourier support also lies on V q . In this section we show that for 
7 G [1/2,1) the local in time solution may be continued for all time, while in the case 7 6 (0, 1/2) the same 
result holds but under the additional assumption that the initial data is small with respect to k. 

Theorem 6.5 (7 > 1/2: Global existence for large data). Let s, ©o and S be as in the statement of 
Theorem \6.3\ If '7 > 1/2, the unique smooth solution of the Cauchy problem (II. 10b — d 1. lib is global in 
time. 

Proof of Theorem lo31 Given the conditions on the initial data and of the source term, and using the prop- 
erties exhibited in Proposition 16.11 the solution constructed in Theorem 16.31 satisfies supp(G) C V q . 
By Lemma 16.21 we have that U is obtained from G by a bounded Fourier multiplier, and hence by the 
Hormander-Mikhlin theorem we have 

||^"||iy s >p < C||0||w a 'P (6.28) 

for all 2 < p < 00, and all s > 0. Therefore the regime 7 > 1/2 becomes "sub-critical" for such solutions, 
since the map i->- U is bounded in Sobolev spaces. Therefore, energy estimates which are similar to 
those used to prove the global regularity of the sub-critical SQG equation (cf. O), may be used with minor 
modifications to prove Theorem l6.5l in the setting 7 > 1/2. 

The case 7 = 1/2 is slightly more delicate. In the case of the critical (SQG) equation, global regularity in 
the critical case has only been established recently (cf. ECOO ). Due to the anisotropic nature of the symbol 
M, it turns out that it is slightly easier to see that the DeGiorgi-inspired proof of Q also takes care of the 
7 = 1/2 case of this theorem. The only fact we must verify is that when G L°°, and supp(0) G V q , 
ther0J7 = M[0] G BMO. Once we prove this, the DeGiorgi iteration scheme of (H, shows that 0(i, •) 
is Holder continuous for t > 0, and one may use energy arguments to conclude the proof of Theorem 16.51 
similarly to iflOl Appendix]. Lastly, to verify that U G BMO, define a new operator M q as the Fourier 
multiplier with symbol M(k)ip q (k), where <p q (k) is a function which is identically 1 on V q , and vanishes 
identically at a fixed distance away from V q . The advantage is that M q is now a periodic pseudo-differential 
operator of order 0, and hence maps L°° to periodic BMO (cf. iMTl ). This is attributed to boundedness of 
M q , inherited from the symbols Mj, and follows from the fact that V q and its collar neighborhood are 
not entirely contained in the curved region of the frequency space where Mj become unbounded. Lastly 
one may verify that when supp(0) C V q , then Af ? [0] = M [0], thereby concluding the proof of the 
theorem. □ 

Theorem 6.6 (7 < 1/2: Global existence for small data). Let 7, s, ©o and S be as in the statement of 
Theorem \6.3\ Ifj < 1/2, there exists a sufficiently small constant e > 0, such that if ||6o||°2 ll©o||/p a + 
|| Qp ||^ 2 1|5|| ^o" oo-jj-s-7 ) ^ £ > where a = 1 — (5/2 — j)/s, then the unique smooth solution of the 
Cauchy problem d 1 - 10b — d 1.11b is global in time. 

Proof of Theorem WEI Since M acts as an operator of order when applied to functions with frequency 
support on V q , from the point of view of energy estimates the situation we are in is exactly the same as for the 
super-critical (SQG) equation. The proof the theorem follows from the same arguments used in ll6l [1411261 
to show the global well-posedness of small solutions to the super-critical (SQG) equation. □ 



'Note that for generic 9 G L°° functions M[Q] does not belong to BMO, but rather to BMO 1 (cf. J9) Section 4]). 
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Appendix A. Smooth solutions to a linear equation 

FOR DATA WITH "THIN" FOURIER SUPPORT 

The goal of this appendix is to prove the existence of smooth solutions to the scalar linear equation 

d t 9 + v -V9 + (-A) 7 <9 = S (A.l) 

9\ t=0 = (A.2) 

where the initial datum 9q, the given divergence-free drift velocity field v, and the external source term S, 
all have frequency support in the same plane V q , i.e. 

supp(<? ),supp(5(V)),supp(?(i, •)) <zV q (A.3) 

for all t G [0, T], for some q G Q, q ^ 0. The main result is: 

Theorem A.l (Existence of solutions with support on a plane). Let 7 G (0, 1), d > 2, ant/ fix s > 
d/2+1-7. Given O G H s {T d ), a divergence-free v e L°°(0, T; F s (T d )), WS G L°°(0, T; iJ s_,y (T d )), 
smc/i f/jaf (IA.3I ) holds, there exists a unique solution 

9 G L oo (0,T;ff s (T ,i ))ni 2 (0,r;F +7 (T d )) 

of the initial value problem (IA. lb — (IA.2b . an J we /zave 

supp(0(V)) CP, 

holds for all t G [0, T]. 

Remark A.2. The main difficulty in proving Theorem I A.l l is in designing an iteration scheme which is both 
suitable for energy estimates, and preserves the feature that in each iteration step the frequency support of 
the approximation lies on V q . In this direction we note that a scheme such that dt9 n+ \ is given in terms 
of v ■ V9 n automatically preserves the frequency support on V q in view of Lemma 16.11 but is not suitable 
for closing the estimates at the level of Sobolev spaces. On the other hand, if we consider dt9^ n+1 ^ to 
depend on v • V9 n+ \, while energy estimates are now clear, it seems difficult to inductively obtain that 

SUpp(6» n+ i) C Vq. 

Proof of Theorem \A. 1\ Since the iteration scheme which is suitable for controlling the frequency support of 
the solution is "loosing" a derivative, we regularize (IA. lb — (IA.2b with hyper-dissipation as 

d t 9 £ + v ■ V9 £ + (-A)W - eA9 £ = S (A.4) 

9 £ (0,-) = 9 (A.5) 

for e G (0, 1], and later pass to the limit e — > in order to obtain a solution of the original system. Since v 
and S are smooth, and v is divergence-free, it follows from our earlier paper O that there exists a unique 
global (or as long as v and S permit) smooth solution 9 s of (1A.4MA.5I ). with 

9 £ G L°°(0, T; H s ) n L 2 (0, T; H s+ ~<) n eL 2 (0, T; H s+1 ). (A.6) 

and the solution is bounded in the above spaces independently of e G (0, 1]. 

The advantage of considering the system (IA.4b — (IA.5b over (IA. lb — (IA.2b . is that for the hyper-dissipative 
system, we can construct a smooth solution (as will be shown below), which has frequency support lying 
in Vq. Since (IA.4b is linear, and we work in the smooth category, i.e. s > d/2 + 1 — 7, the unique 
smooth solution of (IA.4MA.5l l satisfies supp(0 e (t, •)) C V q for all t G [0, T). We note already that the 
e-independent bounds in the regularity class (IA.6b . will allow us to pass to a limit 9 £ — > 9, as e — > 0, and 
this limiting function will automatically satisfy supp(6 l ) C V q , since the latter space is closed and discrete. 

We now proceed to construct a solution 9 s of (IA.4MA.5I ). which has the desired frequency support 
property. We consider the following iterative scheme: the first iterate is given by the solution of 

d t 9\ + (-Ay 9{ - eA9\ = S {A.l) 

0?(O,-) = 0o (A.8) 
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for all t G [0, T), while the n + I s * iterate is given as the unique solution of 

d t 9l +1 +v-V9 e n + (-A)X+i - £A^ +1 = 5 (A.9) 

e e n+1 (o,-) = e (A. io) 

for all n > 1. We note that the solutions of (IA.7MA.8I ) and (1A.9MA.10I ) respectively, may be written 
explicitly using the Duhamel formula 

M -A)+(-Ap)tn , /"* ( e (-A)+(-A)T)(t-T). 



^( t ) = e M-Am-A)')% + / e ^-A W -ArKt-T; 5 ( T)dT (A11) 

JO 

0* (t) = e (-(-A)+(-A)^)t 0() + /" e ( e (-A)+(-A)T)(t-T) (y . (u ^ )(r) + s{t)) ^ (A 12) 

JO 

Since the operator exp((e(— A) + (— A) 7 )i) is given explicitly by the Fourier multiplier with non-zero 
symbol exp ((e|fe| + |fc| 7 )i), this operator does not alter the frequency support of the function it acts on. 
Therefore, it follows directly from Proposition 16.11 our assumptions on the frequency support on 9q an d S, 
that supp(#f (t, •)) C V q for all t G [0, T). We proceed inductively and note that if supp(f^j) C V q , then by 
our assumption on the frequency support of v and Proposition 16. 1 1 we also have supp(t> 9%) C V q . Hence, 

we obtain, as in the case n = 0, that supp(#^ +1 (i, •)) C V q for all t G [0, T), concluding the proof of the 
induction step. This proves that the frequency support of all the iterates 9 e n lies on V q . 

It is left to prove that the sequence {#^}n>i converges to a function 9 s which lies in the smoothness class 
( IA.6I ). Note that there is no cancellation of the highest order term in the nonlinearity: f A s (v ■ \79^ l )A s 9f l+1 . 
However, since at least for now e > infixed, we may use the full smoothing power of the Laplacian. First, 
note that from (1A.4I) it follows that for any t G (0, T] we have 

K x {t) := S up||A s 0?(r,O||i 2 + /V s+7 #i(r, Oilier + e f ^^{r, .)||| 3 dr 
[0,t] J o Jo 

< ||A%||| 2 + f || A s ^S(r, OHiadr. (A.13) 
Jo 

Hence, there exists a time 7\ G (0, T] such that TZ\(Ti) < 2||A S 6> ||| 2 , e.g., any T x such that 

l|A s #oll ?2 

T X < » (A.14) 

I II L oo (0,T;_ff s -T) 

is sufficient. We proceed inductively, and assume that 

K n (t) :=sup||A s ^(r,-)||| 2 + f \\A s+ ^9 n (r, -)\\ 2 L2 dr + e f \\A s+1 9 n (r, -)\\ 2 L2 dr (A.15) 

[0,t] J0 JO 

is such that lZ n (Ti) < 2||A s #o||£2- We now show that if 7\ is chosen appropriately, in terms of e, #o> v an d 
S 1 , we have K n+ i(Ti) < 2||A s #o|||2 too. From (IA.9I ), the fact that V • v = 0, integration by parts, and 
s>d/2 + l — 7>d/2 which makes an algebra, we obtain 

lj t W ASd n+l\\h + ||A^ +1 ||i 2 +,||A s+1 ^ +1 ||i 2 

< ||A>0| i2 ||A s+1 ^ +1 || i2 + ||A s+ ^ +1 || L2 ||A s -^|| L2 

< ^\\\'v\\b\\^Jh + f HA S+1 ^ + illi 2 + ^l|A s+7 ^ + ill! 2 + ^l|A s -^||| 2 (A.16) 
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from which it follows that 



ftn+iCTx) < ||A^o||| 2 + \ [ L \\A s v(T,-)\\ 2 L2 \\A s ef l (T,-)\\l 2 dT+ [ 'llA-^r, Oilier 

£ Jo Jo 

^ II^- S ^o|Il 2 +^1 f~ll U ll£^(O,T;# s )l|A S 0o|||2 + ||'S'|||oo(o i T;ir s -T)^ ■ (A. 17) 



Tl ~~ ZIUI2 \\\sa ||2 , TjTcP (A " 18) 



Hence, it follows from (1A.13I ) that if we let 

e\\V0o\\h 

then we have lZ n+ \(Ti) < 2\\A s 9q\\ 2 l2 - Since the choice of T\ cf. (1A.18I) is independent, of n, it is clear that 
the inductive argument may be carried through, and hence H n (Ti) < 2||A s 6>o||| 2 independently of n > 1. 
To pass to a limit in n, we consider the difference of two iterates, and note that 

9M+1 ~ + ("A) 7 (^ +1 - 91) - eA(9* n+1 -9D + V V(9 £ n - 9^) = (A.19) 

(^+i-^)(O,0=O (A.20) 

for all n > 2. Similarly to (IA. 16b . it follows from dA. 19b that 

K n (t) := sup ||A S (^ +1 - ^)(r, -)||£ a + f ||A s +^(0 n+1 - <£)(r, Ol^dr 
[0,t] «/ 

+ e [ t \\A s+ \9 n+1 -9f l )(T,-)\\l 2 dr 
Jo 

<- /V^(r,-)ll| 2 ||A^(r, OHl^r (A.21) 
£ .In 



and therefore 



R-n{Ti) < ^-||w|li°o( ,T;H s )^-i( T i) (A.22) 
for all n > 2. Thus, due to our choice of T\ cf. (IA. 18b . we have that 

Tl < — (A.23) 

4 \\ V \\l°°{0,T;H s ) 

and hence TZ n (Ti) < H(T\)/2, which implies that the sequence {9 e n \ n >\ is not just bounded, but also a 
contraction in 

L°°(0, Ti; H s ) n L 2 (0, Tr,H s+r ) n eL 2 (0, Ti; F s+1 ). (A.24) 

Hence there exists a limiting function # e in the category dA.241 ). In addition, since for every n > 1 we have 
supp(^) C "P g , and the set V q is closed (and even discrete), we automatically obtain that supp(6> e ) C V q . 

To show that 9 s may be continued in the category (1A.24I ) up to time T, we note that ||A s ^ e (Ti)|| 2 2 < 
2||A s #o|| 2 2 , and hence repeating the above argument with initial condition 9 £ {T\), we obtain a solution 

9 £ e L°°(0,Ti +T 2 ;H s )nL 2 (0,T 1 +T 2 ; H s+ ~f) n eL 2 (0,Ti +T 2 ,H S+1 ), where 

2e]|A s o ]|?2 



OIL, 112 IIAsfl l|2 _i_ <r|| C||2 

°ll c; llL°°(0,T;if 3 )ll-' 1 - "OWL 2 "•" fc ll°llL oo (0,T;_H" s -^) 

which is such that ||A s £ (Ti + T 2 )|| 2 2 < 4||0 O ||| 2 . Since the series 

^ 2 fe e||A^ |li2 



fc> H (; lll,°°(0,T;_ff'»)ll ii U 0\\ L 2 "I" fc IPIl£oo( 0)T ;fl-»-7) 



(A.26) 



diverges for every fixed e > 0, the above argument may be terminated after finitely many steps, concluding 
the construction of the solution 9 £ in the category (1A.6I) . 
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In order to conclude the proof of the lemma we need to pass to a limit as e — > 0. By construction we 
have that 9 £ is uniformly bounded, with respect to e, in L°°(0, T;H s )t~) L 2 (0, T; H s+ ^), and from (IA.4b we 
obtain that dtO 6 is uniformly bounded, with respect to e, in L 2 (0, T; H s ~ 2+1 ). Thus, by the Aubin-Lions 
compactness lemma (see for instance [25]]), we obtain a weak limit G L°°(0, T; H s ) n L 2 (0, T; H s+ ^), 
so that the convergence # e — >• 6* is strong in L 2 (0, T;H S ). Since the evolution is linear and s is large 
enough, it follows that this limiting function is the unique solution of (IA. lb which lies in L°°(0, T; H s ) n 
L 2 (0, T; H s+1 ). Lastly, since for every e > we have supp(# e ) C "P g , and since V q is closed, we obtain 
that the limiting function also has the desired support property, i.e. supp(#) C V q , which concludes the 
proof of the theorem. □ 
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